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The Dyson equation proposed for planar temporal Wilson loops in the context of supersymmetric 
gauge theories is critically analysed thereby exhibiting its ingredients and approximations involved. 
We reveal its limitations and identify its range of applicability in non-supersymmetric gauge theories. 
In particular, we show that this equation is applicable only to strongly asymmetric planar Wilson 
loops (consisting of a long and a short pair of loop segments) and as a consequence the Wilsonian 
potential can be extracted only up to intermediate distances. By this equation the Wilson loop is 
exclusively determined by the gluon propagator. We solve the Dyson equation in Coulomb gauge for 
the temporal Wilson loop with the instantaneous part of the gluon propagator and for the spatial 
Wilson loop with the static gluon propagator obtained in the Hamiltonian approach to continuum 
Yang-Mills theory and on the lattice. In both cases we find a linearly rising color potential. 
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I. INTRODUCTION 



The Wilson loop is a quantity of central interest in Yang-Mills theory. The temporal Wilson loop is related to the 
static potential between infinitely heavy color sources and represents the order parameter of confinement: An area 
law in the temporal Wilson loop corresponds to a linearly rising potential. Furthermore the potential V(L) extracted 
from a rectangular temporal Wilson loop of spatial extension L is the ground state energy of a pair of an infinitely 
heavy quark and antiquark at separation L. As discussed in detail in Ref. [l| this "Wilsonian" potential V(L) is not 
the fundamental color confining potential Vc(L), but rather the residual potential that survives color screening by 
vacuum polarization effects, i.e. V(L) contrary to Vc(L) contains the back reaction of the Yang-Mills vacuum to the 
£N) ' presence of the infinitely heavy color sources. Obviously, the screening of color charges will lower the energy so that 
V(L) < Vc(L). Therefore, a linearly rising fundamental potential Vc(L) does not necessarily imply an area law in 
the corresponding Wilson loop. The fundamental color potential Vc(L) can be explicitly isolated in Coulomb gauge, 
Q\ I Ref. and it was shown in Ref. [|| that the so-called "Coulomb string tension" ac extracted from Vc{L) is indeed 
C^) . an upper bound to the Wilsonian string tension a. 

It is still a big challenge to understand the color confinement mechanism in the continuum Yang-Mills theory, which 
• i-H , requires to prove the area law for the Wilson loop. In recent years several approaches to continuum Yang- Mills theory 
have been intensively pursued. Among others there are the Dyson- Schwinger approach in Landau gauge, Ref. [if, 
and more recently also in Coulomb gauge, Ref. and a variational solution of the Yang-Mills Schrodinger equation 
in Coulomb gauge. Ref. @, 0, H, H, Glj. These approaches have given support for the Gribov-Zwanziger confinement 
scenario, Ref. What is, however, missing in these approaches is the explicit non-perturbative evaluation of 

the Wilson loop, showing the emergence of the area law. In the continuum theory the calculation of the Wilson loop 
is rendered complicated due to the path ordering. In the context of supersymmetric Yang-Mills theory a Dyson type 
of integral equation was proposed for the temporal Wilson loop, which sums all planar ladder (or rainbow) diagrams, 
Ref. [13J. This equation takes care of the path ordering, at least in an approximate fashion, and was recently also 
applied to the temporal Wilson loop in ordinary (non-supersymmetric) Yang-Mills theory, Ref. [3]. Here the only 
input into this Dyson equation is the gluon propagator, which is gauge dependent. In Ref. [l4| the gluon propagator 
was found by solving the Dyson-Schwinger equations of Yang-Mills theory in Landau gauge in the rainbow-ladder 
approximation. In the present paper we will critically review the Dyson equation for the Wilson loop, work out its 
ingredients and shortcomings, and apply it to Yang-Mills theory in Coulomb gauge. We study both the temporal and 
spatial Wilson loop. 

The organization of the paper is as follows: In the next section we give a short heuristic derivation of the Dyson 
equation for the Wilson loop proposed in Ref. [131 ] and critically analyse the ingredients, the approximations involved 
and exhibit its range of validity. In Sect. Ill we solve this equation for the temporal Wilson loop in Coulomb gauge 
assuming an instantaneous, temporal gluon propagator. In Sect. IV the Dyson equation for the Wilson loop with 
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arbitrary gluon propagators is converted to a one-dimensional Schrodinger equation following Ref. [l3j . In Sect. V we 
use this Schrodinger equation to calculate the spatial Wilson loop in the Hamiltonian approach to Yang-Mills theory 
in Coulomb gauge. Finally our conclusions are given in Sect. VI. 



II. THE DYSON EQUATION FOR THE WILSON LOOP 



Below we briefly sketch the derivation of the Dyson equation for the Wilson loop thereby exhibiting the approximations 
involved and working out its limitations. 
Consider the Wilson loop integral 



(W(C)) = 

defined as expectation value via a functional 
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Here C is a closed loop in Euclidean space, d r is the dimension of the group representation and V denotes path 
ordering along this loop. Furthermore A^ = A^T a is the algebra valued gauge field (referred to as gluon field) with 
T a being the anti-hermitean generators of the gauge group and S[A] is the action of the underlying gauge theory. In 
leading order perturbation theory only the quadratic part of the action is kept 



S[A] = \J d d xd d yA^(x)D^(x,y)A u (y) 



(3) 



with Dftvi^x, y) being the bare gluon propagator and d being the number of space-time dimensions, and path ordering 
can be ignored, yielding for the Wilson loop 



W(C) = exp 



-C 2 I{C) 



= 1 - "-C 2 I(C) 



where 



1(C) 



dx^J dy u D tlv (x 1 y) 
c c 



(4) 
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and 



T a T a = -l (d) C 2 (6) 

is the quadratic Casimir operator. For the gauge group SU(N) we have C 2 = ^iv 1 • Obviously, the perturbative 
gluon propagator D(k) ~ p- cannot give rise to an area law. 

Consider now a temporal planar trapezoidal Wilson loop W — W(S, T; L) with two parallel temporal sides of length 
L and T, separated by a spatial distance L, see Fig. [T] This loop consists of four straight paths. Accordingly, the loop 
integral ([5]) receives 16 contributions. In the limit S,T 3> L the dominant contribution comes from the two temporal 
paths (of length T and S, respectively). Keeping only these two paths and ignoring the contributions where both 
integrals in (|5|) run along the same temporal path, this loop integral reduces to 

S T 

1(C) = -2 J dx J dy Q D 00 ((z - Vof + L 2 ) ■ (7) 


The diagram corresponding to this contribution to the Wilson loop is shown in Fig. [2a). Since the gluon propagator 
(exact or perturbative) drops off for large distances one may argue that for S,T ^> L the contribution to the loop 
integral 1(C) (0 with gluon lines connecting the spatial paths as shown in Fig. [D(b) are subleading compared to the 
contribution from the temporal paths shown in Fig. [2Ja). Then in the standard way one may resum all the ladder 
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FIG. 1: Trapezoidal temporal Wilson loop. 
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FIG. 2: (a) Processes involved and (b) processes ignored in the Dyson equation {gj) illustrated in Fig. 



diagrams with gluon lines connecting the two temporal paths shown in Fig. [3]Ja) by a Dyson equation with kernel 
given by the gluon propagator. This yields the integral equation originally proposed in the context of supersymmetric 
gauge theories, Ref. [ljj], which is illustrated in Fig.[3jb) and given by 

S T 

W{S, T; L) = 1 + g 2 C 2 J ds J dtD ((x{s) - s(t)) 2 ) W{s, t; L) , (8) 
o o 

where 

D ((x{s) - x{t)) 2 ) =x-(s)D IMV (x(s),x(t))x+(t) (9) 

and x^j (s) denotes a parametrization of the two temporal paths of the Wilson loop (see the Appendix for an explicit 
realization of these paths). Obviously, the summation of the ladder diagrams by the Dyson equation remains valid 
when the full gluon propagator is used and does not lead to any double counting. 

The above given sketch of the derivation of the Dyson equation clearly exhibits the limitations of this equation: 

1. The Dyson equation (jSJ) can be only applied to strongly asymmetric loops consisting of two opposite long (above 
temporal) and two opposite short (above spatial) paths. Otherwise it does not make sense to include one pair 
of paths (Fig. Ufa)) in the contour integration |5j) and partially sum their ladder diagrams while neglecting the 
contour integrals of the other pair of paths (Fig. [^b)) and their ladders. It is, however, not necessary that the 
Wilson loop is a temporal one. Eq. {8) can be equally well applied to strongly asymmetric spatial loops. Due to 
the limitation to asymmetric loops the static potential extracted from the solution of Eq. (jSJ) will be accessible 
only for spatial distances L<T. In particular, the limit L — > oo is not accessible by this equation. 

2. The solutions of the Dyson equation (jHJ) obviously satisfy the boundary conditions 



W{S, T = 0; L) = 1 , W{S = 0, T; L) = 1 



(10) 
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FIG. 3: Graphical illustration to (a) the Dyson series for the Wilson loop and (b) the Dyson equation 0). 




L L 

T = T — S 

(a) (6) 

FIG. 4: (a) Triangle and (b) rectangular shaped loops representing the Wilson loop W(S,T = 0; L) and W(S, S; L), see text. 

However, for T — or S — the trapezoidal loop degenerates to a triangle shaped loop, shown in Fig. 0] (a) and 
the Wilson loop of the triangle shaped contour is certainly not one. The wrong boundary value (|10p comes with 
no surprise since the limit T = or S = contradicts the assumption S,T 3> L inherent in Eq. ([5]) . The wrong 
boundary value (|10[) has consequences for the Wilsonian potential extracted from the rectangular Wilson loop 
W(S = T, T; L) shown in Fig.[4jb). Assume the solution of the Dyson equation ([8]) yields asymptotically indeed 
an area law as expected for Yang-Mills theory. Due to the boundary value (TTOj) . W(S — T,T;L) misses the 
contribution from half of its enclosed area (i.e. from the triangle shown in Fig. Ufa)). If the remaining triangle 
has to account for the full rectangle the string tension has to be doubled. Consequently we expect from the 
solution of the Dyson equation ([8]) the double of the true string tension. 



3. In the Dyson equation ([5]) the only remnant from the representation of the gauge group is the quadratic Casimir 

Therefore, the Wilson loop obtained as a solution of Eq. ([5]) will show strict Casimir scaling. It is known, 
however, that Casimir scaling occurs only in the intermediate distance regime, Ref. [l5j]. At large distances 
the string tension survives only for the group representations with odd A— ality 1 while the charges of the even 
A— ality representations are screened. Thus again we find that the Dyson equation ([5]) can be appropriate only 
up to some intermediate distances and is bound to fail for large (spatial) distances. 

4. The Wilson loop, per se, is gauge invariant while the right hand side of the Dyson equation ([5]) is gauge dependent 



1 /V-ality characterizes the representation of the center in an irreducible representation of the group. The center of SU(N) is Z(N), 
consisting in the fundamental representation of the N group elements z n = exp (i^j^) l(iv)i n = 0,1,2, . . . , N — 1, Multiplication of a 
group element g of an irreducible representation with 7V-ality k by the center element z n implies the multiplication of g by the factor 
exp(^fc). 
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through the gluon propagator. In fact, the gluon propagator is non-vanishing only after the gauge has been 
fixed. Thus, Eq. © will yield different results for the Wilson loop in different gauges. 

5. Finally, let us also mention that (except for the temporal Wilson loop in Coulomb gauge, see below) the right 
hand side of the Dyson equation J5]) is not renormalization group invariant, while the Wilson loop is. 

From the above analysis it is clear that the Dyson equation ([5]) can be applied only to strongly asymmetric loops 
with two different length scales and that the smaller length scale has to be restricted from above to intermediate 
distances. The limitations of the Dyson equation ([8]) revealed above will manisfcst themselves in the applications of 
this equation to be given below. 



III. THE TEMPORAL WILSON LOOP IN COULOMB GAUGE 



Since the Dyson equation |[H} for the Wilson loop is gauge dependent via the gluon propagator it is preferable to use 
so-called "physical" gauges, for which the gauge constraint and Gauss' law can be explicitly resolved, leaving a gauge 
fixed theory of the physical degrees of freedom. Such a gauge is the Coulomb gauge 

d-A = 0. (11) 

In addition this gauge has the advantage that the quantity gA is renormalization group invariant, which is not the 
case for covariant gauges like Landau gauge. The Dyson equation ([5]) for the temporal Wilson loop contains only the 
temporal component of the gluon propagator 

g 2 D% b (x,y) = {gA* (x)gA b (y)) ■ (12) 

Therefore in Coulomb gauge the Dyson equation ([8} is renormalization group invariant and thus independent of the 
regularization and renormalization scheme. Let us therefore investigate the Dyson equation in Coulomb gauge. In 
Coulomb gauge the temporal gluon propagator has the structure [l[ 

g 2 D a b (x, y) = -6 ab V c (\x - y\)6(x° - y°) + P ab (x, y) . (13) 

Here Vc(|a; — y\) is the so-called non-Abelian Coulomb potential, which is long-range and presumably responsible 
for color confinement. It describes anti-screening of color charges, while P(x,y), also being long-range, describes 
ordinary screening. The latter is responsible for the breaking of the color flux string between the external quarks 
when dynamical quarks are present. Both terms have been calculated in perturbation theory [l|, I17L Il8l Il9j |. The 
non-Abelian Coulomb potential Vci x — y\) was also non-perturbatively calculated both in the Hamiltonian approach 
to continuum Yang-Mills theory 0, 9| and on the lattice [2(1, HJ, H3| ■ At small distance it behaves like the ordinary 
Coulomb potential ~ 1/r, while it rises linearly at large distance ~ <7c T with a coefficient ac, referred to as Coulomb 
string tension, which is somewhat larger than the Wilsonian string tension a extracted from the (temporal) Wilson 
loop (cr c ~ 1-5(7 Ref. H3). 

Since we are mainly interested in the infrared properties of the Wilson loop let us keep here the instantaneous 
(confinement relevant) part of the gluon propagator (|13|) only. This forces us to consider rectangular Wilson loops 

W(T; L) :— W{T, T; L) , (14) 

for which the Dyson equation ([5]) reduces to 

T 

W(T; L) = 1 - C 2 V C (L) [ dtW(t; L) . (15) 
This equation can be converted into the differential equation 



n 



-^W(T-L) = -C 2 V C (L)W{T,L) (16) 

with the boundary condition 

W(T = 0;L) = 1, (17) 
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which is, contrary to the boundary condition of the general case (|T0|) . indeed the correct boundary condition (the area 
enclosed by the loop vanishes for T = 0). The solution of Eqs. (flBj) . (jTTJ) or equivalently of (TT5)) is given by 

W(T; L) = exp {-C 2 V C (L)T) . (18) 

Since the non-Abelian Coulomb potential rises linearly at large distances we have correctly obtained an area law, 
however, with the Wilsonian string tension a replaced by the Coulomb string tension ac > tr. This is expected due 
to the neglect of the non-instantaneous part P(x,y) of the gluon propagator (fT5|) , which screens the non-Abelian 
charge of the static quarks represented by the temporal Wilson lines. Within the instantaneous approximation used 
for the gluon propagator (fT5|) (g 2 Doo(x) — > — Vc(x)S(x )) Eq. (|T5)) is the correct result. It is clear why in the present 
case the Dyson equation © yields the correct Wilson loop: The processes neglected in the Dyson equation © (see 
Eq. J2])(b)) do not exist for an instantaneous gluon propagator. 



IV. EXTRACTING THE STATIC POTENTIAL 

The charm of the Dyson equation © for the Wilson loop is its simplicity. It can be reduced to a one-dimensional 
Schrodinger equation and from the corresponding ground state energy the static potential can be extracted [l3[ . 
Differentiation of Eq. © with respect to S and T yields the differential equation 

92W dSdT L) = d2 ° 2 D + {S - T)2 ) W{S ' T; L) ' (19) 

which together with the boundary condition (fl~0|) is equivalent to the integral equation ((8j) but easier to solve. By 
introducing the variables 

_ S-T 5 + T d 2 1 f d 2 d 2 \ 

L ' L ' dSdT L 2 \dR 2 dr 2 ) 1 J 

this equation is separable [l3[. With the notation W(R, r) :— W(S, T; L) its solution can be expressed as 

W(R, r)=J2 ( c n ex P (VnR/2) + c- exp (-Cl n R/2)) , (21) 

n 

where the (p n (x) satisfy the 1-dimensional Schrodinger equation 



Vn(r) = -^Vn(r) (22) 



with the potential 

U(r) = -g 2 C 2 L 2 D(L 2 (l + r 2 )) (23) 

and the constants have to be chosen such that W satisfies the boundary conditions (JTTJJ) , which in the variables (|2"0)) 
read W(R, r = R) = W(R, r = —R) = 1 for all R. Since the potential is symmetric U(—r) = U(r) the eigenfunctions 
i Pn{i')i n = 0, 1,2, . . . have definite parity ip n (—r) = (—)"</?« (r) and the boundery condition implies that = for 
n— odd. By the symmetry of (f2"Tj) the 0„ can be chosen to be positive definite without loss of generality. 
We are interested in the Wilson loop W(S, T; L) = W(R, r) for L <C S — T — » oo, i.e. r = 0, R oo, which is related 
to the static potential by 

V(L) = -lim hnW(T,T;L). (24) 

1 — >oo 1 

In this limit and for f2„ > the Wilson loop is dominated by the first term in (|2"Tj) with the largest il n for which 
ip n (0) ^ 0, which will usually correspond to the ground state n = of (|2"2")l 

W(T,T;L)\ T>>L = W(R -> oo,r = 0) « Vo (Q)e n ° R / 2 . (25) 
We thus obtain for the static potential (|24|) 

V {L) = -^± + const. (26) 



Thus, to find the static Wilsonian potential we have to solve the Schrodinger equation (|22[) for its ground state as a 
function of L. As we have discussed above the validity of the Dyson-Schwinger equation ([5]) and thus of the expression 
(|2"6")) for the potential is restricted to not too large L. The precise range of validity will depend on the scale set by the 
gluon propagator which defines the potential U(r; L) given in (|23[) . This will be seen in the results. 
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V. SPATIAL WILSON LOOP IN D = 3 + 1 COULOMB GAUGE 



At zero temperature, due to 0(4) symmetry, the Wilsonian potential can be extracted from either the spatial or 
temporal Wilson loops. In the Hamiltonian approach, the spatial loop is, however, more easily accessible than the 
temporal loop, once the vacuum wave functional has been determined. This is because the temporal Wilson loop 
requires also the time-evolution of the vacuum wave functional. As is clear from the discussion given in Sect. 2 the 
Dyson equation ([8]) can be equally well applied to spatial Wilson loops provided the loops are chosen asymmetrically 
with pairwise long and short paths. We are interested in the spatial Wilson loop in the Hamiltonian approach to 
Yang-Mills theory in Coulomb gauge @, H, @| ■ In this case the expectation value is still given by Eq. J2J where the 
functional integration runs now over the transversal spatial components of the gauge field only and the action S [A] is 
defined by 

e~ s W = j(A)\^[A}\ 2 . (27) 

Here J (A) is the Faddeev-Popov determinant and ip[A] is the vacuum wave functional. In the approach of Ref. [tJ 
the variational ansatz for the wave function ip[A] is such that the action functional defined by (f27|) has the form 



S[A] = J d i xd i yAt(x)u J (x,y)Atiy) (28) 

so that the static (spatial) gluon propagator is given by 

D$(x,y) = {At{x)A){y)) = 5 ab ^t i3 {x)uj- l {x,y) . (29) 

Here uj(x,y) is a variational kernel, which is found by minimizing the vacuum energy density, and tij[x) is the 
projector to transversal gauge fields. The result of the variational calculation is confirmed by a lattice calculation 
of the static gluon propagator [23| . see Fig. [5J In the infrared region, which is responsible for the string tension, 
the variational result almost perfectly matches the lattice data. There are only small deviations in the intermediate 
momentum regime (see Fig.[5j, which shows some dependence on the single remaining undetermined renormalization 
parameter [24[ . Possibly, these deviations result partly from the neglect of the three-gluon vertex. One could perhaps 
fine-tune the remaining renormalization parameter to get a better agreement with the lattice data in the intermediate 
momentum regime. However, this region is irrelevant for the asymptotic behaviour of the Wilson loop. The lattice 
data for the kernel w{k) can be nicely fitted by Gribov's formula [ll| 



M 4 
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FIG. 6: The potential U(r; L) i|23D for several distances L. 



where M is a mass scale which, in principle, is determined by the string tension. The fit to the lattice data yields [231 ] 

M fn 880 MeV « 2y/a . (31) 

We will use the static gluon propagator defined by Eqs. (|29p. I|30p. (f3"Tj) in the Dyson equation ([5]) (i.e. in the potential 
(|23|) of the Schrodinger equation (|22[) ) to calculate the spatial Wilson loop in the Hamiltonian approach in Coulomb 
gauge. _ 

When the expression (|30[) is inserted for the (static) propagator D(k) = (2w(fc)) _1 the resulting integrals are UV- 
divergent (see Eq. (|A9[) in the Appendix). The reason is that Gribov's formula does not include the anomalous 
dimension of the gluon propagator. The logarithmic momentum dependence induced by the anomalous dimensions 
is difficult to see on the lattice. In the Hamiltonian approach to the continuum theory p, [7|, which focuses on the 
infrared physics, the anomalous dimension escapes due to the particular variational ansatz chosen for the vacuum 
wave functional, which leaves out the three-gluon vertex and thus the gluon loop contribution to the gluon self-energy. 
The latter is given in this approach by the ghost loop , which dominates in the infrared. The anomalous dimension 
of the static gluon propagator was studied in Ref . [25] . Including the anomalous dimension the UV-term in Gribov's 
formula has to be modified to 



«(*) = y*» (l + aln-±) 2 +^ (32) 

where the anomalous dimension (of the gluon propagator) is given by 7 — 25]. The parameter a depends on 
M. Since the only effect of the anomalous dimension is to make the integrals (|A9[) UV-convergent we will choose 
a small value of a — 0.10. The results do not depend on the precise value of the parameter a. Let us stress that 
the IR-behaviour of the gluon propagator is independent of the anomalous dimension 7. Fig. 7 shows the potential 
U(r; L) (|2"3"|) calculated from the static propagator (|2"9"|) in Coulomb gauge with uj(k) given by Eq. (|3"2")) for various 
spatial distances L. U(r;L) has the form of a double well centered at r — 0. The dip in the potential, necessary 
for the formation of a bound state, flattens as L increases and vanishes for L > 0.5 fm. However, the bound state 
disappears already for L » 0.35 fm. The Wilson potential can then no longer be extracted from the "ground state 
energy" Qo(-L). This limits the use of the Schrodinger equation (f2"2")l to rather small distances. 

The Wilsonian potential V(L) obtained from the ground state eigenvalue flo(L) via Eq. (|26|) is shown in Fig. [7] as 
a function of the spatial distance L. The potential behaves like the ordinary Coulomb potential at small distances 
(below 0.15 fm) and rises linearly in an intermediate distance regime between 0.15 and 0.35 fm. At larger distances 
L the Schrodinger equation (|2"2")l has no bound state and thus Eq. pi?)) ceases to be applicable. The linear rise in the 
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FIG. 7: Left panel: The full static quark potential V(L) obtained from the full propagator J 32)) and the perturbative potential 
V P crt{L) obtained from the perturbative propagator \33\i . Right panel: The full potential minus its perturbative part. 



potential is clearly exhibited when the perturbative potential V pcrt (L), calculated with the perturbative gluon energy 

WpcrtW (l + aln^-^ , (33) 

is subtracted from the full potential, as shown in Fig. [7J Except for the very small distance regime, which is subject to 
numerical inaccuracies, the potential void of its perturbative part can be nicely fitted by a linear function (see Fig. [7J 

V(L)-V peTt (L)=c + aL, (34) 

from which we extract the string tension 

a « (600 Me V) 2 . (35) 

The string tension is about twice its input value of a = (440 MeV) 2 used to fix the scale. This is in accord with our 
remark given in Sect. 2. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the Dyso n eq uation for the Wilson loop proposed in Ref. [IH in the context of supersymmetric gauge 
theory and applied in Ref. [14| to ordinary Yang-Mills theory in Landau gauge. We have examined the approximations 
involved and exhibited its range of applicability. We have shown that this equation can be applied only to strongly 
asymmetric Wilson loops with the smaller extension of the loop limited from above. This equation was then applied to 
the temporal Wilson loop in Coulomb gauge, keeping only the instantaneous part of the temporal gluon propagator. 
In this case the resulting Wilsonian potential coincides with the non-Abelian Coulomb potential. Finally, the Dyson 
equation was used to calculate the spatial Wilson loop in the Hamiltonian approach to Yang-Mills theory in Coulomb 
gauge. Although from the solution of this equation the Wilsonian potential can be extracted only up to intermediate 
distances, a strictly linearly rising potential was found when the perturbative potential was subtracted. The emergence 
of the area law in the Wilson loop found from the static Coulomb gauge propagator is consistent with the perimeter 
law for the 't Hooft loop found with the wave functional of the variational approach in [24| • 

An exact Dyson-Schwinger equation for the Wilson loop is given by the so-called loop equation [26j]. Like all Dyson- 
Schwinger equations this equation is not closed but is part of a tower of integral equations, which need truncations 
and approximations to arrive at a closed manageable set of equations. Such studies should eventually be carried out. 
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APPENDIX A: THE SCHRODINGER POTENTIAL 



The potential (|23[) in the 1-dimensional Schrodinger equation is, up to a constant prefactor, given by the quantity (J9J), 
which is the gluon propagator in coordinate space contracted with the two "temporal" paths x (s). Below we will 
calculate this quantity for transversal gluon propagators in arbitrary dimensions d. 

In Euclidean coordinate space a (translationally invariant) transversal gluon propagator can be expressed by the 
Fourier integral 

(x) = J (8^ - %X) e%kx D (k 2 ) (Al) 

where fc M = k^/Vk? and D (fc 2 ) is a function of fc 2 = k^k^ only. Obviously D^ v (x) has the tensor structure 

D llv {x)=8 llu I{x)-I lxv {x) (A2) 

where 

d d k 



= / 7— l ^ kxD ( fc2 ) (A3) 



and I(x) = Ip^ix). Since x^ is the only external vector, on which I^ v {x) depends, the following tensor decomposition 
holds 

Inv(x) = S/xv^ix) + Xfj,X v I 2 (x) , Xfj, = —f= . (A4) 

vr 

Contracting this equation with 5^ and Xp_x v we find 

dl 1 (x)+l 2 {x) = I(x) (A5) 
I 1 {x)+I 2 (x) = I(x) , (A6) 

where I(x) = x^I tlu {x)x u . Solving the last two Eqs. for J 1 ' 2 (a;) yields 

(d-l)/ 1 ^) = I(x)-I(x) 

(d-l)I 2 (x) = dl(x) - I[x) . (A7) 
Inserting these relations into Eq. (|A4[) we find for the transversal propagator in coordinate space 

(d-l)D llv {x)=5 flu [{d-2)I(x)+l{x)] - x„x u [dl{x) - I(x)\ . (A8) 

The remaining integrals 



I{x 2 )j J (2n) d v ' \{k-x) 

are worked out in the standard fashion using spherical coordinates in fc-space and putting the d-axis of fc-space parallel 
to x^. The integrals over the first d— 2 angles are trivial yielding the volume of the unit sphere Sd-2 in d— 1 dimensions 

(d-l)n— 

= ~ 7TJTV- ■ A10 

r(^) 
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The integrals over the last angle yield with z = x ■ k 

1 



dz 



. izkx 



= 2 



sin kx 



dzz 2 e lzkx = 2 



kx 

sin kx 4 



kx k 2 x 2 



cos kx 



sin kx 
kx 



(All) 



Inserting these results into (|A9j) we obtain 

oo 

I{x) = 2C d { dkk a ~ L D{tf) 



d-iyKfi2\ sinter 



I{x) = 2C d I dkk d ^D{k 2 ) 



kx 

sin kx 2 



kx k 2 x 2 



cos kx 



sin kx 
kx 



where 



a, 



(d- l)7rV 

(2^) d r(^±i) 



(A12) 



(A13) 



To work out the remaining one-dimensional integrals requires the explicit form of the scalar function D (fc 2 ) . For 
definiteness we consider a planar temporal Wilson loop in the 0-1-plane and choose the parametrization of the temporal 
pieces as 



a£(«)= [s,±-,0,- 



For these paths we have 



and 



so that 



and 



x±(s) = (1,0,0,...)=:^, 

x„ :=x+(S)-x-(T)=L(r,l,0,0,...) 
x 2 = L 2 {l + r 2 ) 



(x-e) = 



l + r 2 

With these results we eventually find for the quantity ([9]) from (|Al2j) 

(d-l)D (x 2 ) = [(d - 2) I (x 2 ) + I (x 2 )] [dl (x 2 ) - I (x 2 )] 

where x 2 is given by Eq. (|A17[) . 



(A14) 

(A15) 
(A16) 
(A17) 

(A18) 
(A19) 
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